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SECTION  1 


INTRODUCTION 

The  development  of  many  types  of  low  density,  high  strength 
construction  materials  has  provided  numerous  possibilities  for  efficient 
structural  design.  Quite  often,  however,  the  full  potential  of  these  sophisti- 
cated material  concepts  cannot  be  exploited,  due  to  analytical  shortcomings. 
Economical  and  convenient  analysis  techniques  are  fast  becoming  the  most 
valuable  tools  of  the  structural  designer.  This  observation  is  especially 
true  in  view  of  current  trends  toward  the  development  of  automated  design 
procedures. 

The  practical  applications  of  sandwich  materials,  in  particular,  have 

outdistanced  analytical  capabilities  to  a considerable  extent.  A large  body 

1-4* 

of  literature  is  concerned  with  the  analysis  of  sandwich  structures, 

much  of  it  devoted  to  the  important  question  of  elastic  stability.  However, 

only  the  most  basic  geometries  have  been  studied  extensively.  The  stability 

of  flat  sandwich  panels  $ ^ and  of  full  cylinders  , for  example,  have 

received  much  attention.  A few  investigators  have  reported  on  the  stability 

15-19 

analysis  of  open  cylindrical  shells  , most  often  within  the  context  of 

a linear  theory. 

The  present  report  describes  the  stability  analysis  of  an  open 
cylindrical  shell  of  a somewhat  more  general  type  than  has  been  considered 
previously.  Numerical  solutions  are  compared  with  numerous  analytical 
and  experimental  values.  Some  investigation  is  also  made  regarding  the 
limits  of  applicability  of  the  linear  buckling  analysis  for  sandwich  shells. 


Numerical  superscripts  indicate  References  listed  at  end  of  report. 
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1.1  PROBLEM  DEFINITION 


The  present  analysis  addresses  the  problem  of  incipient  general 
buckling  of  a cylindrical  sandwich  panel  loaded  by  axial  compression  and 
edgewise  shear  forces.  Each  boundary  of  the  shell  is  considered  to  be 
simply  supported,  and  so  constrained  as  to  prevent  shear  within  the  core 
parallel  to  the  panel  edge.  The  component  layers  are  confined  to  linear 
elastic  response  and  small  displacements,  and  are  taken  to  be  rigid  in  the 
direction  normal  to  the  surface  of  the  panel. 

The  stiffnesses  of  the  core  in  extension  and  bending  are  assumed  to 
be  sufficiently  small  that  the  core  strain  energy  is  due  exclusively  to 
transverse  shear  deformations.  The  transverse  shear  properties  of  the 
core  may  be  different,  as  is  the  case  in  many  practical  situations. 

The  sandwich  face  sheets  are  considered  as  thin  composite  shells. 

That  is,  each  face  consists  of  a number  of  orthotropic  layers,  having  different 
constitutive  properties  and  arbitrary  orientations  with  respect  to  the  struc- 
tural coordinate  axes  of  the  panel.  Since  the  flexural  rigidity  of  the  faces 
about  their  own  axes  is  included,  the  analysis  may  also  be  applied  to  the 
buckling  problem  of  a laminated  or  isotropic  shell.  The  face  sheets  are 
assumed  to  be  thin  in  comparison  to  the  radius  of  curvature  of  the  panel, 
and  to  deform  according  to  the  Love-Kirchhoff  assumptions. 

1.2  METHOD  OF  SOLUTION 

The  mathematical  model  developed  in  this  report  is  based  upon  the 
principle  of  minimum  potential  energy.  The  total  potential  energy  of  the 
panel,  consisting  of  internal  strain  energy  and  the  potentials  of  applied 
forces,  is  formulated  in  terms  of  the  displacements  within  the  sandwich 
core.  The  introduction  of  kinematically  admissible  assumed  displacement 
modes  into  the  potential  energy  yields  a quadratic  form  in  certain  discrete 
displacement  parameters.  The  discrete  counterpart  of  the  minimum  potential 
energy  principle,  which  requires  that  the  gradient  of  this  quadratic  form 
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vanish,  produces  a set  of  algebraic  equations  of  equilibrium,  in  the 
familiar  form  of  a generalized  eigenvalue  problem.  The  critical  load  of 
least  magnitude  for  the  continuous  structure  is  approximated  by  the  smallest 
eigenvalue  of  the  discrete  system  of  equations.  The  corresponding  eigen- 
vectors are  approximations  to  the  buckled  displacement  shapes  which  are 
the  eigenfunctions  of  the  continuous  problem. 

The  displacement  formulation  used  in  the  present  analysis  is  a direct 
extension  of  the  sandwich  plate  theory  presented  in  Reference  20.  This 
method  has  been  employed  in  the  buckling  analysis  of  flat  sandwich  panels^ * 
with  good  results. 

1.  3 COMMENTS  ON  THE  LINEAR  THEORY  OF 
STABILITY  OF  SHELLS 

The  classical  stability  theories  of  beams  and  plates  are  founded  upon 
certain  assumptions  which  are  easily  visualized  as  applicable  to  a wide 
class  of  practical  problems.  The  consideration  of  shell  stability  by  a linear 
theory  is  limited,  however,  by  two  particularly  troublesome  effects.  The 
curvature  of  the  shell  surface  couples  transverse  deflections  with  the  in- 
surface displacements.  Consequently,  significant  changes  in  the  surface 
geometry  may  occur  prior  to  the  onset  of  instability,  and  indeed  the  buckling 
phenomenon  itself  is  somewhat  less  clearly  defined  than  for  the  case  of 

plates  or  beams.  Secondly,  the  curvilinear  geometry  tends  to  magnify  the 

1 8 

effects  of  the  slightest  deviations  from  a smooth  initial  configuration. 
Typically,  the  influence  of  such  imperfections  is  most  pronounced  for  shells 
having  severe  curvature,  and  often  the  prediction  of  buckling  by  the  linear 
theory  is  extremely  inaccurate.  When  displacement  coupling  or  irregularities 
are  known  to  be  significant,  the  use  of  a theory  which  considers  these 
effects  is  appropriate. 

The  linear  theory  of  shell  stability,  though  restrictive  by  nature, 
does  provide  an  economical  method  of  analysis  for  problems  in  which  the 
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effects  of  curvature  are  relatively  small.  For  example,  the  shallow  cylindri- 
cal shell  geometry  can  often  be  analyzed  by  a linear  theory  with  good  results, 
since  the  deformations  can  closely  resemble  those  of  a flat  plate. 

The  sensitivity  of  shell  structures  to  initial  geometry  suggests  that 
careful  consideration  be  given  to  the  analytical  representation  of  the  boundary 
conditions.  The  solution  for  simply- supported  edges  is  often  considered  as 
a bounding  case  for  panels  which  are  supported  in  some  manner  along  each 
edge;  however,  the  presence  of  free  edges  or  weak  supports  may  have 
considerable  effect  upon  the  buckling  behavior  of  curved  panels. 

The  limitations  mentioned  in  the  linear  theory  of  stability  of  shells 
should  by  no  means  be  regarded  as  an  indication  that  the  theory  is  without 
value.  Rather,  the  complex  nature  of  the  shell  stability  problem  requires 
that  the  assumptions  inherent  in  the  linearization  process  be  applied  in  a 
judicious  manner,  and  only  when  appropriate,  if  meaningful  results  are  to 
be  obtained. 
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SECTION  2 


DISPLACEMENT  FORMULATION  OF  THE  POTENTIAL  ENERGY 

In  the  formulation  of  the  linear  buckling  problem,  it  is  assumed  that 
two  pertinent  equilibrium  configurations  can  be  identified.  The  initial 
equilibrium  position  is  exactly  that  predicted  by  the  classical  linear  theory 
of  elasticity,  while  the  second  configuration  corresponds  to  some  "adjacent" 
state  which  is  obtained  from  the  first  by  virtue  of  some  small  perturbation. 
The  phenomenon  of  buckling  is  viewed  as  a process  of  sudden  transition 
between  these  two  states  of  equilibrium.  By  adopting  the  first  position  of 
equilibrium  as  a reference  state,  and  considering  only  those  deformations 
which  must  occur  during  the  excursion  toward  the  adjacent  state,  the 
necessary  conditions  for  the  existence  of  the  additional  state  of  equilibrium 
are  defined. 

With  this  visualization  of  the  problem,  it  is  often  appropriate  to  select 
as  a reference  state  that  configuration  for  which  the  prebuckling  displace- 
ments have  already  occurred,  and  the  critical  stresses  are  present,  in 

some  specified  proportions.  In  other  words,  the  problem  can  be  viewed  as 

22 

one  of  initial  stress.  The  critical  load  magnitudes  are  then  taken  to  be 
that  combination  of  the  initial  stresses  for  which  some  additional  equilibrium 
position  can  be  shown  to  exist.  Accordingly,  the  displacement  field  which  is 
considered  represents  that  transition  which  occurs  between  the  two 
equilibrium  states;  that  is,  the  displacements  due  specifically  to  Ruckling. 

In  the  linear  formulation  of  the  problem,  it  is  further  assumed  that  the 
equilibrium  state  existing  before  the  time  of  buckling  is  characterized  by 
displacements  which  introduce  only  negligible  changes  in  the  geometry  of 
the  component. 
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2.  1 GEOMETRY  AND  DISPLACEMENT 

The  cylindrical  panel  geometry  is  shown  in  Figure  1.  The  midsurfaces 

of  the  individual  sandwich  faces  are  selected  as  reference  surfaces,  having 

the  local  coordinates  x , y , z . In  terms  of  the  global  coordinates  r,  0,  x, 

s s s 


(2.1) 


r - R ; s = 1,  2,  c. 


A global  coordinate  z = r - R is  also  employed,  where  R is  the  radius  of 

o o 

curvature  of  the  panel  midsurface. 

The  face  sheets  are  taken  to  be  thin  layers  which  obey  the  Love- 
Kirchhoff  assumptions.  Noting  that  the  Lame'  parameters  of  the  shell  are 
A = 1,  A = r,  the  displacement  at  the  generic  point  within  face  f is  cal- 

l u 

culated  from: 


ZfWf , x 

( Vf  1 
zf  f,y  • 57  1 

f = 1,  2. 


(2.2) 


Since  only  general  instability  is  to  be  considered,  the  variation  in  the  trans- 
verse displacement  through  the  thickness  of  the  panel  is  neglected.  That  is, 

W!  = w2  = w* 

For  the  purpose  of  describing  only  those  displacements  which  are 
due  to  the  buckling  deformation  of  the  panel,  it  is  assumed  that  displace- 
ments which  occur  prior  to  the  onset  of  instability  are  small  and  introduce 
only  negligible  changes  in  the  undeformed  geometry.  The  uncoupling  of 
these  two  sets  of  displacements  is  accomplished  by  postulating  that  there 
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Figure  1 

Panel  Geometry  and  Coordinate  System 


exist  certain  "neutral"  surfaces  within  the  core  layer,  throughout  which  no 
insurface  displacements  are  executed  during  the  buckling  deflection.  In 
particular,  it  is  assumed  that 


(z  - e ) 0 
x 

(z  - e ) $ 

y 


(2.3) 


where  0 and  $ are  functions  of  the  surface  coordinates  x , y , and  e , e 

c c x y 

are  constant  parameters.  Figure  2 shows  the  geometry  of  the  assumed 
displacement  state  (Equation  2.  3)  at  a section  parallel  to  the  straight  edges 
of  the  panel.  The  displacement  u^  indicated  in  the  figure  is  a prebuckling 
displacement  which  is  automatically  excluded  from  the  analysis  by  virtue 
of  Equation  2.  3.  Notice  that  the  functions  0 and  >P  can  be  interpreted  as 
rotations  of  the  core  layer  which  are  executed  about  points  lying  within 
the  cylindrical  surfaces  z = e^  and  z = e^,  respectively. 

The  displacement  fields  in  the  faces  and  core  are  related  by  conditions 
of  continuity  at  the  bond  lines.  Letting  the  total  thickness  of  the  panel  be 

d = t,  + t,  + t , 

12c 

the  equations  of  continuity  can  be  stated  as 


(U  ,V  ,W  ) 
c c c 


= <U2’V2'W2> 


Z = 2 tZ 


(2.4) 


(U  V ,W  ) 
c c c 


= (Uj.Vj.Wj) 


Z = - 2 + S 
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Defining 


and  requiring  that 


cf  = (-D 

Pf  =1  -‘f 

a = P,  + c e 
fx  f ix 

a = P +c  e ; f = 1,  2, 

fy  f f y 


Rf  >:>-L  , 

2 


(2.5) 


the  continuity  conditions  (Equations  2.  4)  can  be  solved  to  yield  the  following 
for  the  midsurface  displacements  of  the  faces: 

“f  = -cf,af*0- 


-cffa£yMtfw-y) 


(2.6) 


w,  = w = w 
f c 


; f = 1,  2. 


The  displacements  of  an  arbitrary  point  within  the  face  sheets  are  then  defined 
by  the  substitution  of  Equation  2.  6 in  Equation  2.  2, 


U£  = -Cf(afx^-2tfw'x)-zfw'J 

2 

Vf  = -C(afy0-ItfW’y)(1+iI) 


!,  W, 

f y 


(2.7) 


W = w 


; f = 1,  2. 


Thus,  the  displacement  components  of  any  point  within  the  sandwich  panel  can 
be  expressed  in  terms  of  the  transverse  deflection  w,  the  core  rotations  0 and 
ll>,  and  the  neutral  surface  parameters  e^  and  e . 

2.  2 STRAIN  AND  STRESS 

The  extensional  strains  and  in- surface  shear  strains  within  the  sandwich 
core  contribute  only  negligible  amounts  of  the  total  elastic  strain  energy  since 
the  corresponding  stresses  are  assumed  to  be  small.  Since  the  core  is  con- 
sidered to  be  incompressible  in  the  direction  normal  to  the  panel,  the  strain 


energy- 


due  to  deformations  in  this  direction  is  also  neglected.  The  remaining 


components  of  strain  in  the  core  are  the  transverse  shear  strains. 


= u + 

W 

XZC  c,  z 

c,  X 

(2.8) 

V 

= V + 

W - — . 

yzc  c,  z 

c,y  r 

Making  use  of  Equation  2.  3,  the  strain-displacement  relations  for  the  core 
can  be  written  as 


y = 0 + w, 

XZC  X 


V=  I T z/R  <ytw’y>' 

7 c 

The  corresponding  shear  stresses  are  calculated  from 

T = G Y 

XZC  XZC  XZC 


T = G Y 
yzc  yzc  yzc 


The  strain -displacement  relations  for  the  face  sheets  are  given  by 

o v 

e . = e , + z * 
xf  xf  f xf 


Cyt  = €yf  + Vyf 


Y f - Y r + r > f = 1,  2. 

xyf  xyf  f xyf 

The  midsurface  strains  and  curvatures  in  Equation  2.11  are  defined  as 


o 

e , = u 

xf  f,  x 


o w 

eyf  = Vf,y  + R. 


Y , = + u , 

xyf  f,  x f,  y 


1 1 


X = -W, 

XI  XX 


X , = -W,  + 

yf  yy 


(2.12) 


X , = -2(w, 
xyf  xy 


where  the  midsurface  face  displacements  uf,  vf  and  w are  defined  in  turn  by 
Equation  2.  6. 

The  generalized  Hooke's  law  equations  for  a laminated  face  sheet  are 
most  conveniently  represented  in  terms  of  the  membrane  stress  resultants  and 
moment  stress  resultants  for  the  entire  face.  The  development  of  these 
equations  is  given  in  Appendix  A . The  resulting  stress- strain  relations  are 
expressed  in  symbolic  form  as 


is!  -ft-fflei 


(2. 13) 


2.  3 POTENTIAL  ENERGY 

With  the  assumptions  of  Section  2.1,  the  reference  configuration  of  the 
panel  is  identified  as  that  state  at  which  the  critical  loads  have  been  applied, 
in  certain  specified  proportions,  but  the  displacements  due  to  buckling  have  not 
yet  occurred.  Although  the  exact  geometry  of  this  reference  configuration  is 
not  known  unless  a classical  solution  is  performed,  its  deviation  from  the 
original  (unstressed)  state  has  been  assumed  to  be  sufficiently  small  as  to  be 
negligible.  Thus  the  applied  loads  assume  the  character  of  initial,  or  residual, 
stresses,  as  noted  previously. 

The  potential  energy  functional  for  an  elastic  body  subjected  to  initial 

stresses  Xa^  and  XT^can  be  expressed  as 
x xy 


^Discussion  and  more  detailed  derivations  of  the  constitutive  equations  for 
orthotropic  laminates  are  widely  available  in  the  literature;  see,  for  example. 
References  24-26. 
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r 


"p  = JIT  2 (°xex  + ayey  + °zez  + TyzYyz  + TxzYxz  + TxyYxy  ) dV 

^m|[oL0)(“4  *<)  (2-14' 

V L 


+ 2 T 


(o) 

xy 


( u,x“ 


, + V,  V,  + w,  w 


x y 


x y 


dV 


The  functions  o ^ ^ and  ^ denote  the  relative  magnitudes  of  the  initially- 
x xy 

applied  stresses,  while  the  parameter  X specifies  the  intensity  of  the  loading. 
For  the  present  problem,  the  region  V includes  both  the  face  sheets  and  the 
core;  making  use  of  this  fact,  and  the  previous  assumptions  concerning  stresses 
and  strains.  Equation  2.  14  is  restated  in  the  form 


2 

TT  = E 


Vyf)dv 


m i ( 


T Y + T Y 
xzc  xzc  yzc  yzc 


) 


dV 


uiJIjj 


f=l  Vy 


N<°>  . N<°> 

xf  2 xyf 

T— W'x  + 2T"  W'xW’y 


dV, 


f 


*f 


(2. 15) 


In  Equation  2.  15  it  is  assumed  that  the  initial  stresses  are  uniformly  distri- 
buted through  the  face  sheet  thicknesses,  and  that  the  squares  of  derivatives 
of  the  insurface  displacements  are  small  in  comparison  with  the  squares  of 
the  transverse  slopes.  By  introducing  the  stress- strain  relations  (Equations 
2.  10  and  2.  13)  into  the  potential  energy,  and  performing  the  integrations 
through  the  panel  thickness,  the  following  is  obtained: 
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dxdy 


It  is  important  to  notice  that,  in  Equation  2.  16,  the  variables  and  limits  of 
integration  refer  to  the  local  surface  coordinates  defined  by  Equation  2.  1. 
That  is. 


x = x 

s 

y = ys  = Rs0, 

b = b = R ©;  s = 1,  2,  c, 
s s 

if  9 = 0 and  0 = 0 and  understood  to  be  the  bounding  straight  edges  of  the  panel. 

The  displacement  formulation  of  the  potential  energy  is  obtained  from 
Equation  2.  16  through  the  substitution  of  the  strain-displacement  expressions. 
Equations  2.6,  2.9,  and  2.  12.  The  resulting  functional  is 


it  = £ f J -ri  ("a  . (a,  0,  - "^t  w,  ) 

p JoJo  llf^fx^’x  2 f ’xx' 


2 ,12 
+ A (a,  0,  - ^t  w ) 

22f  fy  y 2 f yy 


+ A66f(afxN  + <V'x  - tfW'xv) 


xy 


+ 2A 1 2f (afx0>  x " 
+ 2Al6f(afx^’y  + 


bf-zxXttj*- y-  ft"',,1 


V'*'  tfw'«y,fa&0’x  ' IV"’**1 


+ 2A26i'“£>y  + V’*'  Y'VV’-y'  iV'-yy1] 


+ 2cf  [BUf(afx^’xW’xx  " Ivi*1 


+ B2a(V'yw'yy' 


+ 2B66£taf>yw-Xy  + a£y*'x’"'xy  ’ lfw’xy  1 
+ B12£<a£y*-yw’xx  + V’yy  ‘ V“'**w'yy> 


4-  R (CL  (/).  W.  + OC-  0 » W f “ t -W,  w,  ) 
+ 16r  fx^y  xx  £y  x xx  f xy  xx 


* B26£laf>yw'yy+  V'xw>yy  ' t£w’xyw'yy> 

+ 2B  16f^afx^'  xW‘  xy  ‘ it£"'’xxW’xy) 

+ 2B26£lV’yw,-y  ‘ it£w’yyw,-y)] 

+ Sf[vwZ'  Zc£  (A22f  + ) lVV  2lfWW*yy) 


+ (2B22f+  if)  W 1'^’yy’2  • 2B22£w-yy 


+ 2cfD22£(V,yw’yy  ' zVNy* 


* 4(B66f+  if)  (V-x-  il£W-xy)2 


15 


r 


+ 4B.  ,,(a,  - -rt  w,  )(a,  0,  - Tt  w,  ) 

66f  fy  x 2 f xy  fx^  y 2 f xy 


+ 8c  D (a  ip,  w,  - \t  w,2  ) 

f 6of  fy  x xy  Z f xy 


- 2cfA12£fa£>*w  - iV™'**’  ■ 2B12£W- 


+ 2Bi2f,a£x^x-  zV^xx’V'y  ' 2tf“’’yyl 


+ 2c,D,_.(a  ip,  w,  - — t w,  w,  ) 
f 1 2f  fy  y xx  2 f yy  xx 


+ 4B16£(a£x^-x  - It£W>x*"V'x  ‘ iV'-xy1 


+ 4c.D,..(a,  Ip,  W,  - — t w,  w,  ) 
f I6f  fyr  x xx  2 f xy  xx' 


• 2cfA26f(a£x^yW'  zV™'  xy>‘  4B26£WW’: 


- 2°f  (A26f  + 2 ^ ) •V’*’'  ‘ zV^'xy1 


+ 2B_.,(a,  0,  - Tt.w,  )(a  Ip,  - \t  W,  ) 

26f  fx^  y 2 f xy  fy  y 2 f yy 


+ 4c  D , (a,  w w,  ) 

f 26f  fy  y 2 f yy  xy 


+ 2 f 3B26f  + 2"^  ) (V*x  ' iV^xyNVS 


+ 4c, D (a  Ip,  w,  - -t  W,  W,  ) 

f 26f  fy  x yy  2 f xy  yy 
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SECTION  3 


DISCRETE  FORMULATION 


In  this  section  the  potential  energy  (Equation  2.17)  is  cast  into  a 
discrete  form  from  which  a numerical  solution  can  be  obtained  for  the 
buckling  problem.  The  displacement  field  is  expressed  as  a linear  combination 
of  known  functions  which  satisfy"  the  necessary  conditions  of  kinematic  admissi- 
bility. The  corresponding  discrete  form  of  the  potential  energy  is  a quadratic 
function  of  a number  of  parameters:  these  are  selected  in  such  a manner  that 
the  energy  expression  is  minimized. 

3.1  DISCRETIZATION  OF  THE  POTENTIAL  ENERGY 

The  discrete  counterpart  of  a functional  F[u(x,  y)]  is  obtained  by  the 

introduction  of  an  assumed-mode  approximation  of  the  form 

N 

u(x,  y)  = £ c.y.(x,  y),  (3. 1) 

j=l  J J 

where  c.;  j = 1,  2,  . . . , N are  arbitrary  constants,  and  X.(x»  y)  are  known 
J J 

functions  satisfying  the  imposed  boundary  conditions  and  appropriate  require- 

27 

ments  of  continuity  and  differentiability.  For  the  displacement  formulation, 

stress  boundary  conditions  are  satisfied  as  stationary  conditions,  while 

22 

displacement  boundary  conditions  are  of  the  imposed  type.  Thus,  in  the 
present  problem,  the  assumed-mode  approximations  to  the  functions  0(x,  y), 
^(x,  y)  and  w(x,  y)  are  required  to  satisfy  the  displacement  conditions 

0(x,  o)  = 0 (x,  b)  = 0 

4£»  (o,  y)  = t!>  (a,  y)  = 0 (3>  Z) 

w(x  , o)  = w(x,  b)  = 0 
w(o,  y)  = w(a,  y)  = 0. 

The  above  boundary  conditions  correspond  to  simply- supported  edges  and 
vanishing  core  rotation  on  each  boundary  in  the  plane  parallel  to  the  panel 


IS 


edges.  A suitable  approximation  is  therefore  given  by  the  functions 

M,  N, 

1 rrrrx  . nrT  y 

0(x,  y)  = E E °mncos  — sm-^ 

m=l  n-1 


M1  N1 

y)  = T,  I ^ 

m=l  n=l 


. mrr  x nrty 

sin cos  — 

mn  a b 


M, 

, 2 2 . nxnx  . nn  y 

w(x,  y)  = E E W sin sin-— 

7 , . mn  a b 

m=l  n=l 

Introduction  of  Equations  3.3  into  the  potential  energy  (Equation  2.17)  permits 
the  direct  evaluation  of  the  integrations  over  the  area,  as  described  in 
Appendix  B. 

By  inspection  of  Equation  2.  17,  it  is  evident  that  terms  of  the  form 

2 2 

a $ will  appear  in  the  energy  expression  after  the  integrations  have  been 
fx  mn 

performed.  By  Equation  2.5,  such  terms  are  of  higher  order  than  the 

desired  quadratic,  due  to  the  dependence  of  and  a upon  the  neutral  surface 

locations  e and  e . Recall  that  the  parameters  e and  e , which  are  as  yet 
x y x y 

undetermined,  are  dependent  upon  the  displacement  modes;  that  is, 

e = e (m,  n) 
x x 

e = e (m,  n) 

y y 

The  discretized  energy  expression  can  be  reduced  to  the  desired  Quadratic 
form  by  the  use  of  the  substitutions 


h = e $ 

mn  xmn  mn 

k = e 'I' 

mn  ymn  mn’ 


(3.4) 


An  example  of  this  substitution  is  given  in  Appendix  B. 

With  these  conventions,  the  potential  energy  is  expressed  in  the  matrix- 

form 


A11 

A12_ 

X1 

*-xi 

Xj  T 

~B11 

B12l 

X1 

at 

12 

A22_ 

k 

X2 

T 

^2 

B22J 

X2 

(3.5) 


19 


The  construction  of  the  matrices  of  Equation  3.5  is  described  in  Appendix  C. 


3.  2 APPLICATION  OF  THE  PRINCIPLE  OF  MINIMUM  POTENTIAL  ENERGY 

The  variational  statement  of  the  principle  of  minimum  potential  energy 
is  expressed  as 

6 tt  = 0 . (3.6) 

P 

The  continuous  displacement  functions  upon  which  tt^  depends  are  in  turn 
specified  uniquely  by  the  conditions  of  Equation  3.6.  For  the  discrete,  linear 
formulation,  the  potential  energy  rr  is  a quadratic  function  of  a number  of 
undetermined  coefficients  (see  Equation  3.1), 

tt  = tt  (c.;  j = 1,  2,  . . . , N). 

P P J 


The  coefficients  are  selected  by  requiring  that  the  gradient  of  this  function 
vanish: 


dTT 

2-  = 0;  j = 1,  2, . . . , N.  (3.7) 

dc . 

J 

For  the  present  problem,  the  discretized  potential  energy  (Equation  3.  5) 
is  to  be  minimized  with  respect  to  the  components  of  the  vectors  JX^J  and  jX^J 

Thus  it  is  required  that 


A11 

A1 2 ~ I 

X1 

1 

1 4-  X 

i 

Pn 

B12 

T 

A„„ 

t 

A12 

22  _ 

2 

1 

LB12 

22 

X- 


(3.8) 


As  noted  in  Appendix  C,  the  matrices  [B^]  and  [B^]  are  nu®-  matrices.  Using 
this  fact.  Equation  3.  8 can  be  restated  as 


rvix.)  + [ajp2i-  0 

[Ai/lX,!  +[A22]|X2|  ♦ *lB22]|X2p  0 
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(3.11) 


Combining  Equations  3.9  and  3.10  yields 
and 

<[A2Z)  - AtB22]  lx2  [• 

Equation  3.12  has  the  form  of  a generalized  eigenvalue  problem,  which  can  be 
solved  for  the  load  intensity,  X,  and  the  vector  { X£  [of  transverse  displacement 
modes.  The  remaining  mode  shape  coefficients  j X^are  then  determined 
directly  from  Equation  3. 11. 


21 


SECTION  4 


NUMERICAL  RESULTS 

The  numerical  solution  of  Equations  3. 11  and  3.  12  is  compared  with  a 
number  of  existing  analytical  results  in  this  section  to  confirm  its  accuracy. 
Comparisons  are  made  with  experimental  data,  and  the  range  of  validity  of 
the  linear  analysis  procedure  is  discussed. 

4. 1 SOME  COMMENTS  ON  THE  NUMERICAL  SOLUTION 

In  performing  the  numerical  solution  of  Equations  3.11  and  3.12,  the 
number  of  terms  to  be  retained  in  the  displacement  functions  (Equation  3.  3) 
is  generally  dependent  upon  both  the  geometry  of  the  panel  and  the  type  of 
loading  to  be  considered.  A panel  having  an  aspect  ratio  which  is  significantly 
different  from  one,  for  example,  may  require  a large  number  of  terms  for  an 
adequate  representation  of  the  displacement  state.  Likewise,  the  presence 
of  a shear  load  often  dictates  the  use  of  a higher-order  displacement  approxi- 
mation, since  the  deflected  shape  of  the  panel  can  be  quite  complex.  For  the 
results  reported  in  this  section,  critical  compression  loads  are  calculated 
using  M^  = N^  = 3 and  M^  = = 6 in  Equation  3.  3;  solutions  for  critical 

shear  loads  or  combined  loads  are  based  upon  the  limits  of  summation 

M = N =6  and  M = N = 8.  The  use  of  a greater  number  of  terms  in  the 
11  2 2 

series  for  the  transverse  deflection  is  motivated  by  the  fact  that  this  is  the 
predominant  mode  of  deformation;  numerical  experiments  verify  this  fact. 
Furthermore,  increasing  the  limits  M^  and  N^  has  a significant  effect  on 
computation  time,  since  the  matrix  [A^],  which  is  to  be  inverted,  is  of  the 
order  4M^N^  (see  Appendix  C). 

4.  2 COMPARISON  WITH  ANALYTICAL  RESULTS 

Numerical  solutions  obtained  by  the  present  method  are  compared  with 
previously  published  analytical  results  in  Tables  1 to  3.  Geometries  considered 
include  sandwich  plates,  laminated  composite  plates,  and  cylindrical  shells  of 
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both  isotropic  and  sandwich  construction. 

Critical  buckling  loads  calculated  for  flat  sandwich  panels  are  compared 
with  previously  reported  results  in  Table  1.  The  first  example  is  a square 
panel  loaded  in  axial  compression.  The  present  analysis  yields  a value  of  the 
critical  load  which  corresponds  very  closely  with  the  nonlinear  finite  element 
result  of  Reference  28.  The  average  experimental  value  obtained  from  four 
tests  reported  in  Reference  32  is  N = -283.  3 lb.  /inch,  indicating  that  both 
results  are  in  error  by  about  nine  percent.  The  result  of  the  present  analysis 
corresponds  exactly  with  that  of  Reference  21  for  the  second  example  of  Table  1. 
This  is  to  be  expected,  since  the  method  of  Reference  21  is  a special  case  of 
the  analysis  described  here.  The  effect  of  a reduced  order  of  approximation 
for  the  core  rotation  functions  (see  paragraph  4.1)  is  illustrated  by  the  last 
two  problems  considered  in  Table  1.  The  solution  for  a pure  shear  load  using 
= 6 and  = 8 with  the  present  method  agrees  within  one  percent 

with  the  eight-term  solution  of  Reference  21;  for  the  final  example,  involving 
combined  compression  and  shear,  the  difference  in  the  two  solutions  is 
insignificant. 

Table  2 summarizes  the  calculated  values  of  critical  loads  for  several 

orthotropic,  laminated  thin  plates.  The  first  two  examples  cited  are  taken 

from  the  nonlinear  analysis  of  Reference  28.  The  comparison  of  the  present 

results  with  those  of  the  nonlinear  analysis  is  very  good.  The  second  panel 

listed  in  Table  2 also  shows  close  agreement  with  the  experimental  value 

N = -21.  7 lb.  /inch  reported  in  Reference  29.  A thin  plate  subjected  to 
xcr 

a pure  shear  loading  is  the  third  example  of  Table  2.  The  calculated  buckling 
load  is  in  agreement  with  the  five-term  Ritz  method  solution  of  Reference  24. 

Solutions  for  the  critical  loads  of  cylindrical  shells  of  isotropic  and 
sandwich  construction  are  given  in  Table  3.  The  first  case  considered  is 
that  of  an  isotropic  shell,  analyzed  by  a linear  method  in  Reference  30.  The 
results  correspond  exactly  in  this  instance.  For  the  second  problem,  a sand- 
wich shell  loaded  in  axial  compression  is  analyzed.  A considerable  improve- 
ment is  observed  over  the  results  of  Reference  19  in  this  case.  An  analysis 
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TABLE  1 

SANDWICH  PLATES  --  COMPARISON  WITH  ANALYITICAL  RESULTS 


1 2.3  4 


Ef 

9.  5 x I06 

10.  6 x 106 

10.  6 x 106 

9.  5 x 106 

vf 

. 30 

. 30 

. 25 

.30 

t 

f 

. 021 

. 025 

. 025 

. 021 

G 

xzc 

19000. 

50000. 

4359. 

19000. 

G 

yzc 

19000. 

20000. 

4359. 

19000. 

t 

c 

. 181 

. 25 

. 25 

. 181 

a 

23.  5 

60. 

10. 

20. 

b 

23.  5 

20. 

40. 

20. 

Loading 

Compression 

Compression 

Shear 

Compres  sion 
and  Shear* 

Reference 

28 

21 

21 

21 

Reference 

Result 

-307. 5 

-1025. 7 

±1344. 3 

N = -411.  5 

X 

N = ±137.  2 
xy 

Present 

Analysis 

-309.4 

-1025.  7 

±1374. 1 

N = -411.  8 

X 

N = ±137. 3 
xy 

*such  that  N = 3N 

x xy 
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TABLE  2 


THIN  LAMINATED  PLATES  --  COMPARISON  WITH 
ANALYTICAL  RESULTS 


1 2 3 


A11 

vO 

2 

<\l 

o 

r—i 

3.  118  x 105 

3.  027  x 106 

A12 

. 850  x 106 

. 177  x 105 

. 091  x 106 

A22 

1.  021  x 106 

3. 118  x 105 

. 303  x 106 

A66 

. 862  x 1 0^ 

. 273  x 105 

. 207  x 106 

D11 

1.  Oil  x 106 

127. 

2.  523  x 103 

°12 

. 842  x 106 

4.  45 

. 076  x 103 

D22 

3 

1.011  x 10 

30.  2 

. 252  x 103 

°66 

. 854  x 103 

6.  88 

. 173  x 103 

t 

. 109 

. 055 

. 100 

a 

11. 

10. 

20. 

b 

9.  75 

10. 

20. 

Loading 

Compression 

Compression 

Shear 

Reference 

28 

28 

24 

Reference 
Re  s ult 

wm 

-19.  2 

±194. 9 

Present 

Analysis 

m 

-19.  1 

±195. 0 

A 
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identical  to  that  of  Reference  19,  bat  neglecting  the  flexural  stiffnesses  of  the 

face  sheets,  is  presented  in  Reference  33.  By  that  method,  a value  of 

N = -760  lb.  /inch  is  obtained  for  the  second  case  of  Table  3.  The  last 
xcr 

example  considered,  a cylindrical  sandwich  shell  loaded  in  shear.,  corresponds 
almost  exactly  with  the  numerical  value  given  in  Reference  16.  However,  an 
experimental  result  of  N = 227  lb.  /inch,  reported  in  Reference  31,  suggests 
that  the  numerical  solution  is  grossly  in  error.  This  discrepancy  is  due  in 
part  to  the  depth  of  the  shell,  as  discussed  in  the  following  sections,  and  to 
difficulties  which  are  typically  encountered  in  experimentally  maintaining  a 
true  simply-supported  boundary  condition  in  the  presence  of  a uniform  shear 
load. 

4.  3 COMPARISON  WITH  EXPERIMENTAL  DATA 


The  present  analysis  has  been  used  to  predict  the  critical  loads  of 

cylindrical  panels  in  axial  compression  for  a number  of  known  experimental 
13  34 

cases.  ’ A comparison  of  the  analytical  and  experimental  results  is 
summarized  in  this  section. 

An  important  comment  is  necessary  concerning  the  effects  of  differences 
between  the  theoretical  and  experimental  boundary  conditions,  and  the  impli- 
cations of  the  present  choice  of  assumed  mode  functions.  In  the  experimental 
studies  reported  in  References  13  and  34,  the  loaded  edges  of  the  specimens 
were  attached  to  metal  strips,  while  the  unloaded  edges  were  supported  only  by 
loose-fitting  wooden  guides  slipped  onto  the  edges.  Such  a method  of  support 
resembles  the  theoretical  simple  support  provided  the  included  angle  of  the 
panel  is  small.  However,  for  very  deep  panels,  it  appears  that  the  experi- 
mental edge  condition  is  more  flexible  than  the  theoretical  support  condition, 
since  the  analysis  boundary  condition  imposes  a rigid  reaction  in  the 
direction  exactly  normal  to  the  panel.  Both  boundary  conditions  are  pictured 
schematically  in  Figure  3.  In  addition  to  this  apparent  discrepancy  in  the 
boundary  conditions,  it  is  important  to  note  that  the  assumed  mode  functions 
(Equation  3.3)  employed  in  the  analysis  are  periodic  with  respect  to  the 
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TABLE  3 


r 

f' 

i 


CYLINDRICAL  SHELLS  --  COMPARISON  WITH 
ANALYTICAL  RESULTS 


Exf 

1 x 107 

1 x 107 

4.  19  x 106 

Eyf 

1 x 107 

1 x 107 

4.  01  x 106 

vxyf 

. 30 

. 30 

. 150 

Gxyf 

. 385  x 107 

. 385  x 107 

. 562  x 106 

*£ 

. 020* 

. 010 

. 021 

G 

xzc 

N.  A.  ** 

10000. 

32000. 

G 

yzc 

N.  A. 

10000. 

18300. 

t 

c 

N.  A. 

. 1875 

. 300 

R 

50. 

50. 

21.94 

© 

34.  4° 

34.  4° 

1 28? 

a 

30. 

30. 

49. 

b 

30. 

30. 

49. 

Loading 

Compression 

Compression 

Shear 

Reference 

30 

19 

16 

Reference 

Result 

-81.  0 

-823. 

±630. 

Present 

Analysis 

-81.  0 

-773. 

±629. 

*Total  thickness  of  isotropic  shell 

**Not  Applicable 
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IDEALIZED  SIMPLE  SUPPORT 


TYPICAL  EXPERIMENTAL  SUPPORT 


Figure  3 

Comparison  of  Theoretical  and  Experimental  Boundary 

Conditions 
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circumferential  dimension  of  the  panel,  so  that  the  trial  space  of  functions 
disallows  the  accurate  representation  of  deformations  which  are  axisymmetric 
in  nature.  As  the  included  angle  of  the  panel  becomes  large,  the  buckling 
deflections  can  be  expected  to  approach  the  axisymmetric  deformations  which 
are  typical  of  a full  cylinder.  For  these  reasons,  the  present  analysis  should 
be  expected  to  yield  accurate  results  only  for  panels  which  art  sufficiently 
shallow  that  the  form  of  the  assumed  modi  appr  .ns  and  the  chosen 

boundary  conditions  are  appropriate. 

The  test  data  of  Reference  14  art  .. 
numerical  results  obtained  by  1 1 • - <•  < 

specimens  has  aluminum  fat  < st  • • . g 


’ corresponding 

■ ip«  rimental 


in  thickness,  bonded  to  an  t nd-t  <i  = 

* 


.125  inch  . The  core  shear  m e .. 
..  35 


;d  . " 3 2 inch 

I knrss  of 
> r.  >>aa<*d  upon 
material, 

1 1 iui  >nu« » that 

produces 

■ nt . 

t ev.it  c t ed  to  be 


average  moduli  for  th<  dens:'.  ■ :L  ■ 

and  therefore  are  not  exact.  F m 

an  increase  in  the  core  modulus  - mi. 

an  increase  in  the  calculated  critn  it  . i hut 

The  lack  of  precise  data  for  th<  cor.  n ■ .. 

critical  to  the  comparison  of  numi-ri.  .i  . • ' ■ r . . g 1 . ads. 

The  comparison  of  analytical  and  expt-riri  • t it  r.  .its  . T able  4 indicates 

reasonably  good  agreement  for  panels  having  small  included  angles.  Relative 

o 

errors  in  the  calculated  values  range  from  4.  5“..  (for  6 14  ) to  more  than 

175%  (for  G = 42.9°).  Furthermore,  the  least  errors,  for  panels  having  large 
included  angles,  are  observed  for  those  panels  having  the  smallest  overall 
surface  dimensions.  This  trend  is  significant  in  that  the  influence  of  the 
experimental  boundary  conditions  is  indicated,  since  the  smaller  panels  are 
better  represented  by  the  theory. 


1 


*Exact  dimensions  for  each  of  the  specimens  are  tabulated  in  Reference  34. 
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TABLE  4 


COMPARISON  OF  ANALYTICAL  VALUES  WITH  EXPERIMENTAL 
DATA  OF  REFERENCE  34. 


Specimen 

No. 

Length 
a (in.  ) 

Radius 

R (in.  ) 
c 

Angle 
© (deg. ) 

N 

xcr 

(exptl.  ) 

N 

xc  r 
(calc.  ) 

55-1 

29.  0 

63.  2 

17.  4 

-165. 

-229. 

55-2 

29.  0 

63.  2 

17.  3 

-183. 

-233  . 

56-1 

28.  9 

59.9 

14.  0 

-220. 

-230. 

56-2 

28.  9 

63.  2 

13.  2 

-200. 

-220. 

58-1 

41.  1 

51.  0 

46.  3 

-473. 

-687. 

58-2 

41.  0 

48.  3 

48.9 

-430. 

-735. 

58-3 

41.  0 

51.  7 

45.  8 

-487. 

-699. 

60-2 

59.9 

77.  4 

42.  9 

-518. 

-1432. 

61-1 

70.  0 

9 2.  9 

42.  7 

-170. 

-419. 

61-2 

59.9 

104.  0 

32.  6 

-199. 

-353. 

61-3 

69.9 

73.  2 

54.  2 

-199. 

-484. 

All  specimens  Aluminum  faces  on  balsa  wood  core. 

E.  = 10.  x I06  psi,  Vr  = • 30,  G 10000.  psi 
f i c 


30 


1 


The  experimental  data  reported  in  Reference  13  have  also  been  studied 
by  the  present  analysis  method.  The  test  specimens  are  constructed  of 
aluminum  face  sheets,  between  . 012  inch  and  . 032  inch  in  thickness,  bonded 
to  corkboard  cores.  Pertinent  geometrical  data  are  given  with  the  comparisons 
of  critical  loads  in  Table  5.  The  relative  errors  in  the  calculated  buckling 
loads  are  seen  to  be  quite  poor  for  this  group  of  specimens. 

One  important  difference  exists  between  the  two  sets  of  test  specimens 
described  above.  For  the  panels  described  in  Reference  34  the  core  shear 
moduli  are  approximately  10,  000  psi  based  on  average  values  taken  from 
Reference  35.  The  moduli  of  the  test  specimens  reported  in  Reference  13  are 
between  320  and  15,  000  psi,  but  fall  primarily  in  the  range  320  - 1000  psi. 

Thus,  the  analysis  procedure  is  more  accurate  when  the  shear  rigidity  of  the 
sandwich  core  is  relatively  large. 

The  relative  errors  observed  for  numerically  determined  buckling  loads 
in  axial  compression  are  plotted  versus  panel  included  angle  in  Figure  4.  This 
data  includes  the  cylindrical  panels  from  References  13  and  34,  and  two  flat 
panels  (see  paragraph  4.1)  cited  from  References  29  and  32.  A number  of 
conclusions  can  be  made  from  this  plot: 

1.  For  panels  having  very  mild  curvature  (included  angles  of  approxi- 
mately 15°  or  less),  the  linear  buckling  analysis  gives  reliable  results 
in  most  instances.  The  error  committed  in  the  linear  analysis  tends 

to  increase  very  rapidly  past  this  limit,  however.  A similar  limita- 
30 

tion  is  reported  for  isotropic  shells,  for  included  angles  larger 
than  about  30°. 

2.  The  error  in  the  numerical  solution  is  the  most  serious  for  sandwich 
panels  having  relatively  weak  cores.  This  suggests  that  the  prebuckling 
deformation  may  be  more  significant  for  weak  core  sandwich,  or  that 
the  adoption  of  a linear  displacement  field  through  the  core  thickness 

is  unduly  restrictive  for  such  problems. 
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TABLE  5 


COMPARISON  OF  ANALYTICAL  VALUES  WITH 
DATA  OF  REFERENCE  13 


Specimen 

No. 

9114-17 

61-1 

61-3 

1155-10 

9114-10 

1463-25 

1155-25 

9114-25 

1155-6 

1463-3 

1155-3 

9114-6 

9114-3 


Length 
a (in.  ) 

29.9 
70.  0 

69.9 
30.  0 
29.  7 

29.4 

29.9 
29.  8 

29.9 

29.4 

29.9 

29.  8 

30.  0 


Radius 

R (in.  ) 
c 

23.  8 

92.9 
73.  2 

24.  4 
26.  2 

22.  9 

23.  2 
27.  2 
60.  5 
23.  6 
27.  7 
58.  6 
26.  0 


Angle 
©(deg.  ) 

29.  2 
43.  0 
54.  4 
28.6 

26.4 
29.  8 

29.8 

25.  2 

18.9 
28.6 
25.  2 
19.  5 

26.9 


N 

xcr 

(exptl.  ) 

-373. 

-170. 

-199. 

-466. 

-267. 

-779. 

-776. 

-479. 

-255. 

-384. 

-289. 

-132. 

-225. 


EXPERIMENTAL 

N 

xcr 
(calc. ) 

-1186. 

-421. 

-486. 

-1467. 

-773. 

-2652. 

-2024. 

-1096. 

-475. 

-1305. 

-942. 

-292. 

-586. 
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ERROR  IN  CALCULATED  BUCKLING  LOAD 


I 

! 


3.  Scatter  observed  in  the  experimental  values  increases  rapidly  with  the 
panel  included  angle.  This  fact  supports  the  previous  comments 
regarding  degeneration  of  the  experimental  boundary  conditions  for 
large  included  angles. 
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SECTION  5 


i 

■ 


SUMMARY  AND  CONCLUSIONS 

A method  of  analysis  has  been  presented  for  the  prediction  of  stability 
of  open  cylindrical  sandwich  shells  subjected  to  axial  compression  and  shear 
loading.  The  mathematical  model  is  obtained  from  the  linearized  principle 
of  minimum  potential  energy.  The  total  potential  energy,  consisting  of 
internal  elastic  strain  energy  and  applied  loads  potentials,  is  formulated  in 
terms  of  displacements  within  the  core  of  the  sandwich.  The  Ritz  method  is 
used  to  transform  the  resulting  functional  into  a quadratic  function  of  a set  of 
discrete  parameters,  which  are  then  selected  in  such  a way  that  the  energy  is 
a minimum.  The  infinite  sequences  of  eigenvalue  and  eigenfunctions  for  the 
continuous  problem  are  approximated  by  the  finite  number  of  eigenvalues  and 
eigenvectors  thus  obtained,  the  lowest  eigenvalue  corresponding  to  the  load 
intensity  for  which  instability  is  first  observed. 

The  analysis  method  described  has  been  compared  with  numerous 
existing  results,  both  theoretical  and  experimental.  Agreement  with  the 
experimental  results  is  shown  to  vary  in  a somewhat  predictable  pattern  as  a 
function  of  the  included  angle  of  the  panel,  and  to  be  influenced  by  both  the 
overall  panel  dimensions  and  the  shear  rigidity  of  the  core.  At  least  some  of 
the  disagreement  between  theory  and  experiment  can  be  attributed  to  the  varia- 
tion in  experimental  edge  conditions,  although  further  study  is  required  in 
order  to  evaluate  the  effect  accurately.  A few  theoretical  values  obtained  for 
shallow  panels  with  relatively  strong  cores  are  encouraging,  showing  relative 
errors  of  five  to  ten  percent  based  on  test  data. 

The  present  study  has  demonstrated  that  while  the  linear  theory  of  sand- 
wich shell  stability  is  restricted  to  a rather  small  class  of  geometric  configura- 
tions, it  is  possible  to  obtain  useful  results  by  respecting  the  limitations  of  the 
theory.  For  deep  shell  geometries,  it  is  evident  that  a more  sophisticated 
analysis,  incorporating  the  effects  of  prebuckling  deformation,  is  required. 
Further  study  will  help  to  define  the  limitations  of  the  linear  formulation  with 
greater  clarity. 
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Thus,  Equation  A2  relates  the  stresses  and  strains,  referred  to  the  structural 
axes  of  the  shell,  in  layer  k of  the  face  sheet. 

By  Equation  2.11,  the  strains  are  separated  into  uniform  and  linearly- 
varying  components. 
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where  z is  the  distance  from  the  midsurface  of  the  face, 
combining  Equations  A2  and  A3  is  written  symbolically  as 
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The  stress  and  moment  resultants  for  the  face  sheet  are  defined  to  be 
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where  the  superscript  k refers  to  the  particular  layer  at  a station  zf  within 
the  integration.  Integration  of  Equation  A4  through  the  thickness  with  respect 
to  weighting  functions  1 and  z^.  yields 
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The  elements  of  the  3x3  symmetric  matrices  A^,  and  D are  defined  by 
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By  Equations  A6  and  A7,  the  membrane  stress  resultants  and  moment  resultants 
are  related  to  the  midsurface  strains  and  curvatures  in  a simple,  concise 
manner. 
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appendix  b 


INTEGRATION  FORMULAS 


The  following  formulas  are  used  in  performing  the  area  integrations 
of  the  discretized  potential  energy  function: 
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Here  i and  j are  integer  indices,  and  6..  is  the  Kronecker  delta. 

6 -l1’1"3  (P2) 

ij'  ( 0 ; i^j. 

The  operator  A is  defined  in  such  a way  that,  for  any  algebraic  function  X, 

°dd 

lJ  ( o ; i + j even. 

The  above  definition  is  understood  to  apply  also  to  the  case  where  the  function  X 
is  indefinite,  as  for  i = j in  the  last  two  integrals  of  Equation  Bl. 

The  mechanics  of  performing  the  integrations  is  best  illustrated  by  a 
simple  example.  Consider  the  first  term  of  the  potential  energy  (Equation2.  17), 
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Substituting  the  assumed  mode  approximation  (Equation  3.3)  gives 
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Next,  the  integrations  are  performed  according  to  Equation  Bl.  Using 
Equations  B2  and  B3,  the  resulting  expression  is 
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Substituting  for  a.  from  Equation  2.  5,  one  obtains 
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It  should  be  noted  that  the  terms  exmn*mn  and  exmn  mn  are  exPresslons  of 

cubic  and  quartic  order,  respectively,  for  reasons  described  in  Paragraph  3.1 

Ilsinc  the  definition  for  h (Equation  3.4),  terms  of  the  desired  quadratic 
6 mn 

order  are  obtained: 
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The  collection  of  the  above  and  similar  expressions  into  a convenient  matrix 
form  is  outlined  in  Appendix  C. 
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STIFFNESS  AND  GEOMETRIC  MATF-.IX  FORMULAS 


Definitions  are  made  in  this  Appendix  for  the  entries  of  the  elastic  and 
geometric  stiffness  matrices  of  Equation  3.5.  The  symbols  5..  and  A.  are 
used  as  defined  in  Appendix  D.  The  following  notational  convention  is  also 
extablished:  The  superscript  T,  when  assigned  to  an  expression  F(i,  j,  k,  1,  ) 

is  taken  to  mean 

j,  k,  1)  = F(k,  1,  i,  j), 

where  i,  j,  k and  1 are  integer  indices. 

The  assumed-mode  coefficients  (Equation  3.3)  are  arranged  in  order 
in  the  vectors 
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The  assembled  vectors  of  unknowns  (Equation  3.5)  are  then  defined  as 
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Note  that  the  dimension  of  |X^|is  4 M^hy  and  that  of  jX^Jis  M^N^. 


The  elastic  stiffness  matrix  is  defined  by 
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where  the  numbering  A*^  refers  to  the  ordering  of  the  vectors  of  assumed- 

1 3 

mode  parameters  within  { X^}  and  { X For  example,  A is  the  submatrix 
obtained  from  the  products  of  $ ^h  in  the  potential  energy  expression. 

kl  kl 

The  individual  entries  A_  of  the  submatrix  A are  given  by 
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where  the  following  definitions  are  used: 
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For  convenience  in  the  numerical  solution  of  the  assembled  matrix 
equations,  the  elastic  stiffness  matrix  is  partitioned  into  blocks  corresponding 
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This  form  of  the  stiffness  matrix  corresponds  to  that  used  in  Equation  3.5. 


Definitions  of  the  entries  of  the  geometric  stiffness  matrix  are  made 
in  a similar  fashion.  However,  due  to  the  assumption  that  the  squares  of  the 
derivatives  of  the  insurface  displacements  are  small  in  comparison  with  the 
derivatives  of  the  transverse  slopes  (see  paragraph  2.3),  the  submatrices 
B and  are  simply  null  matrices.  The  nonzero  entries  of  B2?  are 
calculated  from 
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where  i and  j are  related  to  m,  n,  r and  s as  noted  previously.  It  should  be 
noted  that  the  expressions  N^^and  N^° ^ refer  only  to  the  relative  magnitudes 
of  the  applied  loads,  since  the  intensity  of  each  load  is  characterized  by  the 
single  parameter  X (Equation  2.  15). 

The  entries  of  B^  may  be  restated  in  terms  of  a single  parameter 
for  each  load  type  by  making  the  definitions 
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